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We report kelvin to millikelvin-temperature measurements of dissipation and frequency shift in 
megahertz-range resonators fabricated from ultra-nanocrystalline diamond. Frequency shift Sf/fo 
and dissipation Q _1 demonstrate temperature dependence in the millikelvin range similar to that 
predicted by the glass model of tunneling two level systems. The logarithmic temperature depen- 
dence of Sf/fo is in good agreement with such models, which include phonon relaxation and phonon 
resonant absorption. Dissipation shows a weak power law, Q _1 oc Ti, followed by saturation at 
low temperature. A comparison of both the scaled frequency shift and dissipation in equivalent mi- 
cromechanical structures made of single-crystal silicon and gallium arsenide indicates universality 
in the dynamical response. 
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I. INTRODUCTION 

Micromechanical resonators are instrumental in the in- 
vestigation of a wide variety of fundamental physics prob- 
lems. These include quantum measurement and quan- 
tum computation [TJ [2], ultra-sensitive force and mass 
detection [3], single spin detection fjQ, gravitational wave 
detection [T] and other fundamental phenomena [5J [5J. 
The relevant features of a resonator are characterized by 
the resonance frequency shift 5f (compared to the low- 
est temperature resonance frequency) and dissipation (in- 
verse quality factor) Reduced dimensions are neces- 
sary for achieving high resonance frequencies. However, 
miniaturization beyond the sub-micron scale leads to a 
dramatic increase in the surface-to-volume ratio, result- 
ing in increased dissipation, limiting device performance. 
To counter this trend it is therefore necessary to avoid ex- 
trinsic mechanisms and minimize intrinsic mechanisms. 

From the low temperature response of single crystal 
silicon resonators at both kilohertz [7] and megahertz [8 
frequencies, it is clear that low-lying energy excitations 
of internal defects or two-level systems (TLS) provide 
the dominant contribution to intrinsic dissipation in mi- 
cromechanical resonators. Even though the general trend 
of the low temperature dependence is easily explained by 
the standard glass model of TLS [3 HO] , additional exper- 
iments in other materials such as silicon [8] and gallium 
arsenide [TT] and detailed theoretical calculations [TH [T3] 
suggest an incomplete understanding of the temperature 
dependence of the quality factor. 

Here, we report a detailed set of low temperature mea- 
surements in a novel material, ultra-nano-crystalline dia- 
mond (UNCD). We find that the resonance frequency 
shift and quality factor of resonators fabricated from 
UNCD, single-crystal silicon, and epitaxially-grown gal- 



* Electronic address: mohanty@bu.edu 



lium arsenide show universality in their temperature de- 
pendence. The universal behavior provides further evi- 
dence that intrinsic dissipation in micromechanical res- 
onators is dominated by crystal defects, whether they are 
(substitutional) crystal impurities or (configurational) 
surface defects caused by the abrupt lattice termination. 
Our measurements indicate that the same TLS processes 
dominate in a range of materials, underlining the need 
for a single material independent theory. 

Polycrystalline diamond is an exciting material for mi- 
cromechanical devices due to its unique combination of 
physical and electrical properties. The UNCD used here 
does not maintain the high thermal conductivity of sin- 
gle crystal and NCD (nanocrystalline diamond), and the 
measured Young's modulus is low compared to single 
crystal and NCD diamond materials. In contrast, UNCD 
is distinguished by the ability to be grown to a wide 
range of thickness without compromising surface fric- 
tion and roughness. NCD the crystal size is generally 
on the same order as the film thickness and grows in 
size with film thickness. UNCD maintains low surface 
roughness through a re-nuclearization process that in- 
duces new crystals to form. The resulting conglomera- 
tion of 5 — 10 nm sized crystals is extremely smooth up 
to considerable thicknesses. In both NCD and UNCD 
electrical properties can be varied over a wide scale by 
controlling the doping concentration. Here, low-doped 
wafers are used, and the material is electrically insulat- 
ing. 

We find that dissipation follows a power law in tem- 
perature, Q~ x cx Ts, independent of resonance frequency 
and material. At low temperatures, below ss 100 mK, 
dissipation saturates. Upon cooling, a shift in the sound 
velocity Sv/vq can be observed. This is equivalent to 
measuring the frequency shift Sf/fo [14] . Below a char- 
acteristic temperature T rs 1-2 K, where resonant- 
absorption processes of TLS are active, we observe a log- 
arithmic temperature dependence with a positive slope. 
Above T<>, TLS relaxation-absorption mechanisms dom- 
inate. The temperature dependence is also logarithmic, 
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however with a negative and steeper slope. Interestingly, 
the slope of log(<5// /o) scales with f$ . We find that 
data from sub-micron-sized mechanical resonators from 
UNCD, single crystal silicon and gallium arsenide het- 
erostructures can be scaled onto a universal curve for 
both Sf/fo and Q' 1 . 
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FIG. 1: (color online) a) Frequency-length -2 relation in- 
cluding a linear fit from which the Young's modulus of 
440 ± 70 GPa is determined. Insert: SEM micrograph 
of the diamond harp structure with three doubly clamped 
beams, b) Displacement vs. force plot for all five fundamen- 
tal resonances, the spring constant is given by the inverse of 
the slope. The arrows mark the onset of non-linearity. Insert: 
22.55 MHz resonance with lorentzian fit of the 11.3 [im beam. 



previously measured in polycrystalline diamond [T7] , 
It is common that thin film diamond contains a large 
number of sp 2 bonds on grain boundaries, which degrade 
the Young's modulus. UNCD is particularly sensitive to 
this as the high density of nano-crystal grains effectively 
increase the area of the surface where the sp 2 bonds 
form. Compressive strain within the film may arise with 
cooling do to the thermal expansion mismatch of poly- 
crystalline diamond and the silicon base. The frequency 
shift caused by tensions is discussed by Postma et 
al. [T5] . Figure [T] d) shows the oscillation amplitude on 
resonance, x(/o) = QFd r /k e ff as a function of the drive 
force, F dr — BLIdr is the actuation voltage) for all five 
beam resonances, where the measured magnetomotive 
response, V em /(cJ ) = gL m f Q -frfr(^o), is proportional 
to the beam amplitude [TS], where B is the applied 
perpendicular magnetic field, m the effective mass of the 
resonator, ujq the resonance frequency and I dr the drive 
current across the beam. £ is an integration factor. The 
effective spring constant k e ff is determined from the 
inverse of the slope of the plot in figure [T] b) and found 
to be on order unity for all structures. For the following 
data, care is taken to drive the resonators only in the 
linear regime. 

III. RESULTS 



II. EXPERIMENTAL SETUP 

The micromechanical resonators (Fig. [l] a) are fab- 
ricated using standard micro-machining procedures (e- 
beam lithography, metalization, reactive ion etching and 
finally an HF wet-etch to release the structures). The 
samples are etched from 340 nm thick Aqua25 films 
obtained from Advanced Diamond Technologies, Inc. 
Mechanical resonances from five doubly-clamped beam 
structures of 350 nm width and lengths varying from 
8.7-19.6 /im are studied. The structures are actuated 
and detected by the magnetomotive technique |15j in 
high vacuum. Resonance frequencies /o and dissipation 
Q^ 1 = A///o are obtained from a lorentzian fit (illus- 
trated for the 22.55 MHz resonance in Fig. [l]b), A/ is 
the full width at half maximum). 

The length dependence of frequency is de- 
picted for the five structures in Fig. [l] c). In- 
clusion of the metallic electrodes on top of 
the diamond modifies the resonance frequency: 
fo = Vy/Wd + E m I m )/(p d A d + p m A m )L- 2 [TS], 
where Ei,Ii,pi, and Ai are the Young's modulus, the 
second moment of inertia, the density and the cross sec- 
tional area respectively of both the diamond and metal. 
L is the beam length and rj a numerical factor. The pre- 
cise fit validates the use of the thin-beam approximation. 
From the slope, we determine the Young's modulus of 
the thin film to be 440 GPa. This is considerably lower 
than expected, when compared to the Young's modulus 



There are a number of dissipation mechanisms that 
are known to contribute to loss in micromechanical res- 
onators. Figure [2] a) depicts the temperature dependence 
for three resonance frequencies; the insert illustrates how 
dissipation varies with frequency at two applied magnetic 
field strengths. Short, doubly-clamped beams often suf- 
fer from clamping losses and dissipation is dominated 
by energy loss into the pads [17]. For this case, dissi- 
pation rapidly grows with decreasing resonator length. 
We observe the opposite behavior, and hence, we can 
rule out clamping loss as a dominant dissipation mech- 
anism. Thermoelastic damping (self-heating) decreases 
with frequency; this mechanism however has been shown 
not to contribute at low temperatures and small struc- 
ture sizes [H]. Internal losses in the electrode also con- 
tribute to dissipation: Q- 1 = [1/(1 + /?)] {QJ 1 + (3Q^ U ) 
with (3 = (EAutAu) I '{E d t d ) — 0.02, where ti is the film 
thickness of the gold and diamond. A conservative es- 
timate for losses in the gold electrode predicts an effect 
smaller than 5%. It can be shown that magnetomotive 
damping (a dissipation mechanism caused by the inter- 
action of the readout curcuit and an externally applied 
magnetic field) [20 is proportional to /q 1 ' 5 - The insert 
of Fig. [2] a) shows good agreement with this prediction for 
applied fields of 6 T and 3 T. At high frequencies, magne- 
tomotive damping no longer contributes significantly as 
indicated by the saturation of the data. For the follow- 
ing dissipation data discussed in this letter, dissipation 
measurements are taken at 3 T, in the magnetic field 
regime where magnetomotive damping can be neglected. 
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FIG. 2: (color online) a) Q 1 — T dependency for three resonance frequencies at 3 T. We measure Q 1 oc T a , where a « |. 
Below 100 mK dissipation saturates. Insert: Dissipation-frequency relation at 6 T and 3 T, including fits for magnetomotive 

damping (Q MMD oc / 2 ) and clamping loss (Q Ci oc / 2 ). Tmc is the mixing chamber temperature, b) Sf/fo — T dependency, 
fo is the resonance frequency at 35 mK. Below and above a characteristic temperature T<> the shift is logarithmic as predicted 
by TLS models. The insert depicts the slopes of df/fo and f ~ x fit for low (up triangles) and high (down triangles) temperature 
regimes, c) Normalized dissipation Q _1 /Q^ for UNCD, Si [8], and GaAs [TT] for frequencies ranging from 5 MHz to 45 MHz. 
The temperature has been scaled to T (see panel b)). The black dotted line represents a slope of ~, alluding to the same 
power law for all three materials. The other slopes represent predictions of various TLS models, d) (<5///o)jv vs. T/T a . fo 
is the frequency measured at 0.1 T . The curves collapse if the ratio s%/si is the same. The dashed lines illustrate four — 
ratios. ^ = —0.5 and |* = — 1 depicts the predictions of the TLS standard tunneling model [5] and the Phillips model |10] 
respectively. 



For the frequency shift data, measurements are taken at 
6 T, as no significant field dependence of this variable is 
observed. 

The presence of TLSs become apparent through the 
temperature dependence of both dissipation, and, in par- 
ticular, resonance frequency, displayed in Fig. [2] a) and 
b) respectively. Dissipation saturates below 100 mK 
with a possible recurrence beginning at temperatures 
above 2-4 K, observable for the lowest frequency. Within 
this temperature range, dissipation follows a power law, 
Q^ 1 oc T a with a m 0.35, and no observed frequency cor- 
relation. At higher field strengths, where magnetomotive 
damping and possibly magnetic field-TLS coupling oc- 
curs, the temperature dependence becomes weaker, while 
the decreasing trend is strongest for lower frequencies 
(data not shown). The lack of a clear frequency depen- 
dence of a measured at 3 T indicates that the measure- 



ments were taken in the parameter space where field cou- 
pling can be neglected. 

The frequency shifts vs. temperature plots in Fig. [2] 
b) show a logarithmic dependence, where at T the slope 
changes sign and increases in magnitude. The logarith- 
mic behavior is in good agreement with glassy TLS mod- 
els, however the slope is expected to remain constant in 
magnitude above T (Phillips model [TO]). The fit results 
are summarized in Table [I] We observe that the magni- 
tude of the slope s, doubles for T > T independent of 
frequency, where si is the slope for T < T and S2 is 
the slope for T > T . The slopes show an inverse de- 
pendence on frequency, also not predicted by standard 
theories. This behavior has been observed by Kleiman 
et al. [7 in previous work on single crystal silicon res- 
onators as well. T shows a rising trend with increasing 
resonance frequency. A precise dependence could not be 
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determined from the available data. 

TABLE I: Fit results for Sf/f - T dependency for both 
temperature regimes (T < T si, T > T ss). The dia- 
mond and silicon resonators show roughly similar slope ratios 
— S2/S1 ~ 2, while for gallium arsenide the magnitude of this 
ratio significantly greater — s%/si ~ 4. The S2 values quoted 
for silicon are rough estimates. T is defined in Fig.[2]b). Er- 
rors on the frequency are on the order of 10 Hz. The error in 
T is considerable due to the limited data points, especially 
above T . This uncertainty is estimated to be ±0.25 K, hence 
the scatter is explained by experimental error. 



Material 


fo 


si 


S2 


_£2 
SI 






[MHz] 


xl0~ 5 


xl0~ 5 




[K] 


UNCD 


5.11 


21.2 


-44.3 


2.09 


1.27 


UNCD 


13.50 


12.0 


-24.8 


2.07 


1.86 


UNCD 


22.55 


7.5 


-19.0 


2.53 


1.61 


UNCD 


30.71 


8.0 


-19.9 


2.48 


1.91 


UNCD 


39.59 


7.2 


-15.3 


2.12 


1.71 


UNCD 


44.85 


7.1 


-15.1 


2.11 


1.92 


Si 


12.03 


1.8 


(-4.6) 


(2.59) 


0.54 


Si 


14.59 


2.0 


(-3.1) 


(1.59) 


0.64 


GaAs 


15.82 


13.5 


-55.7 


4.13 


1.29 



IV. DISCUSSION 

At low temperatures, TLS are confined to a nearly de- 
generate ground state, modeled by an asymmetric double 
well potential with asymmetry A and tunnel splitting en- 
ergies Ao, typically on the order of 1 K. Glass models gen- 
erally assume a constant density of states, independent 
of A and A , given by P(A, A )dAdA = ^dAdA . 
This model gives an expression to determine the number 
of defects from the coefficient of the logarithmic depen- 
dency of the frequency shift with temperature, however, 
it is assumed that this slope is frequency (length) in- 
dependent, and hence does not apply, ft may still be 
useful to use this as a first approximation, where one ob- 
tains on the order of 400 defects for the shortest to 3200 
defects for the longest beam resonators. The Phillips 
TLS model is a more constrained model with a fixed 
splitting energy and a Gaussian asymmetry distribution 
with width Ai, while the density of states takes the 
form P(A)dA = Aexp(-A 2 /2A 2 1 )dA. For both mod- 
els the logarithmic dependence of the frequency shift is 
reproduced, with different coefficients. The Soft poten- 
tial Model (SPM) assumes a quartic potential and con- 
stant density of states for each mode, based on the con- 
cept of localized low frequency sound waves. The SPM 
is successful in explaining low temperature thermal con- 
ductivity and predicts a Q" 1 a Ti dissipation law with 
low temperature saturation |13j . Standard single crystal 
models predict non-logarithmic temperature dependence 
of the frequency shift and non-power-law temperature de- 
pendence of the dissipation, in contrast to the data pre- 
sented here. These theoretical results are summarized 
in [5]. Considering the large surface to volume ratio of 



UNCD one expects a high density of TLSs. Theories have 
been developed to take interacting TLSs into account [9] , 
these too however fail to explain the observed behavior. 

The temperature dependence for both frequency shift 
and dissipation scales to a universal curve for all the 
three materials studied in this temperature range. For 
scaling, the measured temperature is referenced to the 
characteristic temperature T defined in Fig. |2]b). All 
dissipation plots are normalized to the dissipation value 
measured at T , Q~ 1 /Q^ and plotted vs. T/T<>. Fig- 
ure [2] c) illustrates the power law Q^ 1 cx for the 
three materials. Low temperature saturation is observed 
in diamond as well as in silicon, however this is not seen 
in the gallium arsenide resonators for the lowest tem- 
peratures measured. Correspondingly the ^ « —4 for 
gallium arsenide, is significantly greater from that mea- 
sured in diamond and silicon. To compare the shift in 
frequency for different temperatures, the data are nor- 
malized to (<5///o)jv, defined in Fig. [2] b) and plotted 
vs. T/T . For these plots fo is defined as the frequency 
at 0.1 T . Although there are not enough temperature 
data points for the silicon resonators to make statements 
with high confidence, it appears that the three materials 
behave qualitatively the same, as above T the slope of 
<5///o not only changes sign but also increases in mag- 
nitude (see Fig. [2]d)). Such observations are consistent 
with data in literature for a wide range of structures and 
frequencies [71 [51 [TP]. 

The large surface-to- volume ratio introduces a signifi- 
cant amount of defects, concealing the presence of a crys- 
tal; hence the micromechanical resonators behave like 
glassy structures. The data sets appear to fall into uni- 
versality classes defined by the S2/S1 ratio; this must still 
be confirmed with future experiments. The universal be- 
havior of these micromechanical devices motivates the 
need for a theory that captures more accurately the ob- 
served relations. A rigorous theoretical framework must 
be developed that encompasses the density of states of 
TLSs, surface effects due to the high surface to volume 
ratio of micro-electromechanical resonators, and possibly 
a non-ohmic thermal bath. It remains to be seen if, with 
a variation of current TLS models, it is possible to the- 
oretically reproduce simultaneously both the dissipation 
power law and the logarithmic frequency shift. 

While dissipation has been studied extensively at lower 
frequencies in single crystal materials such as silicon and 
gallium arsenide, we present first measurements down to 
millikclvin temperatures of micromechanical resonators 
fabricated from a nanocrystalline material. Understand- 
ing phonon-dynamics (phonon coupling, role of disorder, 
crystal structure) in diamond is fundamentally impor- 
tant to novel condensed matter phenomena such as su- 
perconductivity in boron doped diamond and the pur- 
suit of macroscopic quantum states in mechanical oscil- 
lators [TJ [21]. A necessary first step is to understand 
phonon dynamics in undoped diamond. What is funda- 
mentally exciting is that such micro-resonators demon- 
strate the same universal glassy behavior as those made 
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of silicon and gallium arsenide. 

In conclusion, we report the presence of TLSs in UNCD 
diamond resonators at megahertz frequencies. Both the 
resonant absorption and relaxation absorption regimes 
are observed with a characteristic temperature around 
1-2 K. A constant power law in dissipation and analo- 
gous logarithmic frequency dependence on temperature 
is found and compared to measurements in silicon and 
gallium arsenide micromechanical resonators. Further 



theoretical work is required to fully explain this universal 
behavior, possibly by taking more involved TLS density 
of states into account. 
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